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SUMMARY 

Representative examples are presented of ap- 
plications and development of advanced Ccmputa- 
tional Fluid Dynamics (CFD ) codes for aerodynamic 
design at the McDonnell Douglas Corporation (MDC) . 
Transonic potential and Euler codes, interactively 
coupled with boundary layer computation, and solu- 
tions of slender-layer Navier-Stokes approximation 
are applied to aircraft wing/body calculations. An 
optimization procedure using evolution theory is 
described in the context of transonic wing design. 
Euler methods are presented for analysis of hyper- 
soni c configurations, and helicopter rotors in 
hover and forward flight. Several of these pro- 
jects have been accepted for access to the Numeri- 
cal Aerodynamic Simulation (NAS) facility at the 
NASA-Ames Research Center. 


1 . INTRODUCTION 

Computational fluid dynamics (CFD) is used as a 
routine design technique in the divisional engi- 
neering technology and project organizations of the 
McDonnell Douglas Corporation (MDC). Several 


groups are also active in research and development 
of CFD methods, both at the divisional companies in 
the context of their primary applications, and in a 
more generic sense at the McDonnell Douglas 
Research Laboratories (MDRL) . Figure 1 shows a 
simplified corporate organization chart. 

In addition to the computing facilities at di- 
visional locations, a Cray X-MP/14 is available as 
a central resource under the management of the 
McDonnell Douglas Aerospace Information Services 
Company. When authorized, U. S. Government large- 
scale computing facilities are used for specific 
tasks; cooperative projects with NASA centers have 
been especially fruitful. 


The following examples describe CFD applica- 
tions to aerodynamic design and optimization, and 
new methods under development for computation of 
aircraft, helicopter, and missile flowfields. 

These projects are supported by the McDonnell 
Douglas Corporation Independent Research and Devel- 
opment Program; for several projects, the coopera- 
tion of the NASA NAS Program and/or of Cray 
Research, Inc., is appreciated. 
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Fig. 1 McDonnell Douglas Aerospace Group components active in CFD. 


^Director-Research , Flight Sciences Department. 
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2. Simulation of Power Effects on Complex Configu- 
rations with an Invisci d/Viscous Interaction Scheme 
Based on the Solution of Euler and Inverse Boundary 
Layer Equations (T. Cebeci and L. T. Chen, Douglas 
Aircraft Company (DAC), Long Beach, CA) 

The objective of this project (fig. 2) is to 
combine transonic flowfield codes based on the 
Euler equations with a three-dimensional inverse 
boundary-layer method. Complex wing/body and aft 
fuselage/nacelle/pylon configurations, including 
power effects, will be computed with different grid 
topologies. The results will be compared with 
available experimental data. The computer 
resources used in this project include the NAS 
Cray-2 . 

Successful computation of flows with both 
shock-induced and trailing-edge separation requires 
the combination of viscous/inviscid interactive 
schemes with the Euler equations which provide a 
realistic simulation of the vorticity distribution 
downstream of strong shocks. In the inverse 
boundary-layer method of Cebeci, et al . ( 1986), the 
edge velocity is treated as an unknown and the in- 
teractive solutions are obtained with the Hilbert 
integral approach. The coupling of the inverse 
method with two-dimensional Euler equations is pre- 
sented by Chen, Li and Chen (1987); Chen and Chen* 
extend this approach to three dimensional transonic 
flows over wings and wing/body conf igurat ions . 

The computational procedure of Chen and Chen* 
also extends the hybrid mapping/numerical grid- 
generation method of Chen, Vassberg and Peavey 
(1985) to a new scheme where the mesh lines are 
wrapped around the wing tip into a C-mesh which 
allows proper modeling of wing-tip vorticity 
permitted by the Euler equations, as shown in 
figure 3. 


Task 1: Develop and validate Euler methods for wing/body 

and aft-fuselage/nacelle/pylon configurations, 
respectively, on NAS computer. 

Task 2: Develop and validate an interactive inverse boundary- 

layer method for coupling with the Euler method 
wing/body and aft-fuselage/nacelle/pylon 
configurations, respectively. 

Tack 3: Validate solutions for both configurations on the 
NAS computer for a wide range of flow conditions 
and power settings, and study the mesh effect on the 
solutions by comparing them with DAC full-potential 
and Euler-correction solutions. 

Task 4: Develop a zonal method for coupling solutions about 

wing/body and aft-fuselage/nacelle/pylon 
configurations. 

Task 5: Compare computed solutions with available test data. 
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Fig. 2 Simulation of power effects on complex configurations with 
an inviscid/viscid interaction scheme based on the solution of 
Euler and inverse boundary-layer equations. 


*Chen, L. T. and Chen, H., "An Interactive Scheme 
for Transonic Wing/Body Flows Based on Euler and 
Inverse Boundary-Layer Equations," Douglas 
Aircraft Company, Long Beach, CA, 1987. 



(a) KGRID (C-H-H mesh) 



(b) PGRII) (C-C-H mesh) 



(c) QGRID (C-C-C mesh) 
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Fig. 3 Comparison of cross-sectional grid distribution for three 
grid-generation methods. 


The three-dimensional inviscid flowfield is 
calculated with Jameson's multi-grid, four-stage 
Runge-Kutta time-stepping scheme, and the inverse 
boundary-layer method of Cebeci et al . (1986) is 
applied with a strip theory approximation. Both 
displacement surface and blowing velocity 
approaches are used in the interactive matching 
procedure described in detail by Chen and Chen.* 

Figure 4 shows an inviscid Euler calculation of 
supersonic flow about a fighter-type wing/fuselage 
configuration. The calculation was obtained with 
the new grid-generation method described by Chen, 

V assberg, and Peavey (1985) and Chen and Chen.* A 
comparison of viscous/inviscid interactive calcula- 
tions with test data is shown in figure 5 for a 
high- aspect -ratio transport wing/body configura- 
tion. Both potential and Euler codes were used for 
the inviscid flowfield, and the interactive cou- 
pling was accomplished by the displacement surface 
approach. The agreement between the calculations 
and experiment is generally good, with some deteri- 
oration toward the wing tip. 

Further improvements of the grid-generation 
method will include implementation of the zonal 
technique for complex configurations. 


3. Transonic Wing Optimization By Evolution Theory 

(R. D. Gregg, Douglas Aircraft Company, Long Beach, 

CA) 

The present comments are based on Gregg and The procedure involved 640 FLO-22 solutions 

Misegades (1987), which contains a more extensive with a total CPU-time of 7.2 hours. This example 

discussion of the impact of increasing computer is the first one of three in this paper where the 

capability on three-dimensional optimization proce- microtasking procedure of Booth and Misegades 
dures and presents detailed practical results. (1986) was used in cooperation with Cray Research 

on their Cray X-MP/48, resulting in this case in 

The importance of computing speed is illustrat- clock-time reduction by a factor of 3.95 to 1.8 
ed in figure 6, where the assumption of an optimi- hours, 
zation process requiring 1000 transonic wing flow- 
field solutions indicates the dramatic effect of 

computational speed: the approximately 250 hours 4. Accurate, Efficient Prediction Method for 

required by an IBM 3081 can be reduced to ten hours Supersonic/Hypersonic Inviscid Flow (A. Verhoff and 
on a single-processor Cray X-MP, and to 2.5 hours p. j. O'Neil, McDonnell Aircraft Company, St. 

when taking full advantage of the multitasking cap- Louis, M0) 

ability of a Cray X-MP/48. 

The new method employs a spatial marching pro- 
A description of the evolution- theory approach cedure to solve the formulation of Euler equations 

is shown in figure 7. Basically, a set of merit described by Verhoff and O'Neil (1984). This for- 

functions related to performance requirements is mulation is written in terms of Riemann-type vari- 

calculated repeatedly for different sets of design ables by use of a local streamline coordinate sys- 

variables that are allowed to vary within a chosen tern shown in figure 9, and therefore models wave 

range. At each step, the range and its center are propagation in a physical sense with no inherent 

adjusted until a reasonable convergence of the Mach-number restrictions. The procedure can be 

merit functions is achieved. extended to three-dimensional flows in a straight- 

forward manner, and it is especially suitable for 
Figure 8 illustrates how the convergence of efficient vector coding, 

merit functions was achieved in the case of an 

aspect ratio of ten, taper ratio of four, 25-degree Steady-state solutions are obtained by an ex- 

swept NACA 0012 wing which was allowed to change plicit time-integration scheme, where the local 

camber in the last 15 per cent of chord. The flight maximum time step is used to increase rate of con- 

condition was M=0 .76 at C =0.55, and the base wing vergence. Spatial derivatives are approximated by 

exhibited excessive shock L drag, buffet, and tip one-sided finite differences to properly model the 

stall. While induced drag increased slightly be- wave propagation. At solid surfaces, the algebraic 

cause of increased departure from an elliptic span- relationship between the flow angles 0 and <f>, corn- 

wise load distribution, significant gains in other bined with the local pressure gradient, yields an 

merit functions were achieved in approximately accurate boundary condition. Neither non-physical 

*Chen , L. T. and Chen, H. , "An Interactive Scheme 
for Transonic Wing/Body Flows Based on Euler and 
Inverse Boundary-Layer Equations," Douglas 
Aircraft Company, Long Beach, CA, 1987. 



Heavy line indicates isovalue = 1.0000 
Dashed line indicates isovalue = 1.3000 
Increments in isovalue = 0.0500 


Fig. 4 Mach number contours on a typical fighter wing/fuselage; 
= 1.3 and a = 4.84°. 

forty iterations, including a 50-per cent decrease 
in total drag. 
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input parameters, artificial dissipation, nor 
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Fig. 6 Aerodynamic wing optimization: the importance of 
computational speed. 


• User selects range of each design 
variable to initiate scheme 


_Ro- 


Ao 

(Center) 


A random sample within this range 
is selected and evaluated 


■Initial Range - 


For fully inviscid supersonic flow, three- 
dimensional flowfield analysis can be reduced to a 
sequence of two-dimensional problems by use of a 
spatial marching procedure based on this formula- 
tion. Flowfields about vehicles with wings, tails, 
and fins can therefore be predicted with relatively 
coarse computational grids. Shocks are automati- 
cally captured; the only constraint requires the 
outer grid boundary to include the perturbed flow- 
field. 

A computer program called SCRAM (Streamline 
Coordinate Riemann Axial Marching) has been devel- 
oped at McDonnell Aircraft using this formulation. 
The program evolved from a research code and within 
one year became a useful design method for routine 
analysis of hypersonic configurations. A detailed 
description of the method is presented by Verhoff 
and O’Neil,* including the following representative 
e xam pi es . 

Figure 10 illustrates the suitability of the 
SCRAM-code for vector processing in the case of a 
simple 20-degree cone at zero angle of attack in 



Fig. 8 Aerodynamic wing optimization for a simple, swept-wing test 
case (in terms of induced, viscous, shock, trim, buffet, and 
total drag). 


• Best value of design variable(s) 
is selected as new center, and 
the radius is modified 

* M # I f I 

4 2 1 3 5 6 7 8 

Best - 1 , worst - 8 

• Procedure is repeated until 
‘convergence’ of all design 
variables is obtained 

Amax old 


Amin old 


Ao 


An , 


Ama: 



Fig. 7 Evolution theory. 
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Fig. 9 Definition of flow angles and streamline coordinates. 


*Verhoff, A. and O’Neil, P. J., ’’Accurate, Efficient Prediction Method 
for Supersoni c/Hypersoni c Inviscid Flow,” McDonnell Aircraft Company, 
St. Louis, MO, 1987. 
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M=2.5 supersonic flow; CPU time reductions by an 
order of magnitude can be achieved for higher grid 
densities relative to scalar operation. Figure 11 
shows the accuracy of the solution for various grid 
densities, as compared with the exact analytical 
solution. 

Figure 12 shows the computational grid for a 
blended wing/body configuration. The grid has 
half- plane distribution of 51 by 25 nodes, with 25 
axial stations, and is based on a nonlinear grid 
generation procedure developed at McDonnell Air- 
craft. For M=6 , computed pressure contours in four 
cross-planes are shown in figure 13, and a compari- 
son of predicted force and moment characteristics 
with test data is shown in figure 14; estimated 
skin-friction drag has been subtracted from the 
test data, yielding excellent agreement with the 
inviscid computation. Each case, including input- 
output processing, required only about five minutes 
on the Cray X-MP/14. 



Grid dimensions in cross-plane 

GP6I-1439-10-R 

Fig. 10 Vector efficiency of SCRAM code on CYBER 205 
computer; M O0 = 2.5, 6 = 20°, and a = 0°. 
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Fig. 11 SCRAM solution for circular cone; Moo = 2.44, a-0°, and 
6 = 20 °. 


Figures 15 and 16 show the grid and the com- 
puted Mach number contours, respectively, for a 
complex conical body at M=10. The excellent grid 
resolution and computational efficiency contributed 
to the achievement of a converged solution at the 
first attempt, again requiring approximately five 
minutes on the Cray X-MP/14. 

McDonnell Aircraft is applying the SCRAM code 
to high-quality prediction of increasingly complex 
flowfields, including realistic fighter and hyper- 
sonic aircraft configurations. 



Station 1 Station 2 Station 3 Station 4 
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Fig. 12 Computational grid for AMI-X blended body; vertical tails 


Surface pressure distributions 



Normalized buttline 


51 X 25x24 grid 



Station 1 Station 2 Station 3 Station 4 
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Fig. 13 SCRAM pressure predictions for AMI-X blended body; 
Moo =6 and 5 = 0°. 
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a Test data 

(less skin friction drag) 





Fig. 14 SCRAM force and moment predictions for AMI-X blended 
body; M oo = 6.0. 
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Fig. IS MARCHG grid generation for SCRAM code. 
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Fig. 16 SCRAM code application to complex hypersonic vehicle; 
Moo = 10 and a=0°. 


5. Integrated Flowfield Analysis Methodology for 
Fighter Inlets (R. R. Cosner, McDonnell Aircraft 
Company, St. Louis, MO) 

A methodology to predict viscous flow over in- 
tegrated fighter forebody- inlet combinations has 
been developed at McDonnell Aircraft Company; ef- 
fects of Mach number, Reynolds number, angle of 
attack and engine mass flow are included in the 
prediction. A detailed description of the method 
is given by Cosner (1986), including the following 
examples . 

The procedure is based on the velocity- 
splitting method which provides a relaxation solu- 
tion to the steady-state Navier-Stokes equations 
[Cosner (1982)]. A multiple-zone mesh achieves the 
geometric flexibility required for representation 
of realistic inlet/forebody configurations. The 
mesh is generated by a three-dimensional implemen- 
tation of Thompson’s method, with wall orthogonal- 
ity enhanced by a three-dimensional extension of 
Sorenson’s technique, as described by Cosner 
(1982) . 

All terms, except time derivatives, were re- 
tained in the Navier-Stokes equations; turbulence 
effects were represented by the Cebeci-Smith al- 
gebraic model. Non-conservative and quasi- 
linearized discretization was used, and type- 
dependent differencing, according to the local flow 
velocity, ensured numerical stability without addi- 
tion of artificial viscosity terms. 
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Separate mesh and flow solutions were executed 
in each of the three zones; the flow solutions are 
coupled by matching boundary conditions between 
adjacent zones during the iterative solution proce- 
dure. The mesh solutions were executed to a maxi- 
mum of 30 iterations using an alternating-direction 
implicit scheme; formal convergence was not found 
to be necessary for achievement of a valid mesh. 

The definition of the three computational zones 
is shown in figure 17. The external flowfield up- 
stream and downstream of the inlet highlight com- 
prises the forebody and centerbody zones, respec- 
tively; the internal zone extends from the inlet 
entrance to the throat or the engine face. Mesh 
parameters, such as node distribution or stretching 
coefficients in each coordinate direction can be 
prescribed entirely independently in each zone. 

The flow solution is coupled across the interface 
by interpolating Dirichlet and Neumann boundary 
condition data; mesh coordinates and dependent var- 
iables are modeled using piecewise bilinear func- 
tions. The flow solution algorithm selects the 
appropriate boundary condition according to the 
type of local flow in the adjoining zones at or 
near the interface. 

The iterative solution procedure uses an 
alternating-direction implicit scheme. The cycle 
updates each zone in turn, starting with the fore- 
body. In the internal zone, an algebraic correc- 
tion to provide mass continuity in the potential 
field significantly enhances the rate of conver- 
gence. Typically, 10 to 30 solution cycles are 
required; a larger number is necessary when the 
free-stream Mach number approaches one from either 
direction or the viscous interactions are strong. 
Also, a modified procedure is used at near-choking 
conditions in the inlet. 



□ A^/A^. —0.71 (100% power) 
A A o /A c =0.53 (flight idle) 



Application of the zonal forebody-inlet code to 
the F/A-18 configuration used 8000 computational 
nodes (16x25x20) in the forebody zone, 4000 nodes 
(8x25x20) in the centerbody zone, and 6000 
(12x25x20) in the inlet; the quas i-cyl indri cal mesh 
topology is shown in figure 18. Pressure distribu- 
tions on the forebody in front of the inlet are 
shown in figure 19 at M=0.8 for three engine mass- 
flow rates from near-maximum to zero flow. The 
mass-flow rate variation was introduced into the 
computation only in the inlet zone, from which the 
perturbation propagated through the zone boundary 
into the forebody zone. As shown in Cosner (1986), 
especially good agreement with wind tunnel tests 
was obtained in the critical region inside the 
lower intake lip. The computations were performed 
on a CDC CYBER 176, typical times being 10 to 15 
minutes. 


\ 



\ 

\ 

\ 
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Fig. 17 Zone layout for forebody/inlet analysis. 


Fig. 19 Effect of engine airflow on F/A-18A forebody pressure; 
Mach = 0.8 and a = 0° (A c = inlet mass flow area, 

A c = inlet capture area, corrected to free-stream 
conditions). 


Figure 20 shows the three-zone mesh interfaces 
for an AV-8B forebody/ inlet configuration; a total 
of 53760 nodes were used. Computations made on a 
Digital Equipment Corporation VAX 1 1/780 computer 
required 20 to 25 hours of wall-clock time. Com- 
puted surface-pressure distributions and comparison 
of the lower inlet lip inboard pressure with a test 
at M=0.9 are shown in figure 21. Figure 22 shows, 
for M=0 .67 , the inlet entrance flowfield from which 
detailed total and static pressure distributions in 
the inlet can be derived for design purposes. 

The zonal procedure for solving the Navier- 
Stokes equations is being extended to integrated 
external compression inlets at all Mach numbers of 
interest, and to include the effects of wing and 
canard surfaces. 


6. Euler Calculations of a Helicopter Rotor in 
Hover and Forward Flight (R. K. Agarwal and J. E. 
Deese, McDonnell Douglas Research Laboratories, St. 
Louis, MO). 

The progressive increase of computing power 
over the last decade has advanced the computation 
of helicopter rotor flowfields from solution of 
transonic small-disturbance equations first to that 
of the full potential equations, and recently, to 


116 









codes based on the Euler equations. An Euler code, 
designated MDROTH , has been developed at the 
McDonnell Douglas Research Laboratories for cal- 
culating the flowfield of a multibladed helicopter 
rotor in hover and forward flight. 

The code solves the three-dimensional Euler 
equations in a rotating coordinate system on body- 
conforming curvilinear grids around the blades. 
Euler equations are recast in absolute flow vari- 
ables so that the absolute flow in the far field is 
uniform, but the relative flow is non-uniform. 
Equations are solved for the absolute flow vari- 
ables by employing Jameson's f i nite-volume explicit 
Runge-Kutta time-stepping scheme. Rotor wake ef- 
fects are modeled in the form of a correction ap- 
plied to the geometric angle-of-attack along the 
blades. This correction is obtained by computing 
the local induced downwash with a free-wake analy- 
sis program. The details of the methodology are 
described by Agarwal and Deese (1986). 

A set of test calculations was performed by 
Agarwal and Deese* to verify the code for a model 
rotor in hover and forward flight at various col- 
lective pitch angles. The model rotor has two 
untwisted, untapered blades of aspect ratio equal 
to six and the NACA 0012 airfoil section. Computa- 
tions were performed on a 97 (chordwise) x 33 
(blade-normal) x 21 (spanwise) mesh. 

The code has been fully vectorized for optimum 
performance on Cray X-MP, and has also been micro- 
tasked for peak performance on the f our-processor 
Cray X-MP/48 [Booth and Misegades (1986)] with a 
speed-up factor of 3.71 . A typical case requires 2 
million 64-bit words of main memory and 7x10 
seconds of CPU time per mesh point for each itera- 
tion. A solution for hover flowfield converged in 
500 to 800 iterations. This project has been ac- 
cepted for access to the NAS Cray-2. 

Figure 23 shows the main features of the flow- 
field of a two-bladed rotor in hover, with the im- 
bedded finite-difference grid and the coordinate 
system. The flowfield is characteri zed by tran- 
sonic shocks, complex vortical wakes, and blade- 
vortex interactions. 

Comparisons of computed pressure distributions 
with experimental data are shown in figure 24 for a 
hovering rotor at a tip Mach number of M =0.52 and 
zero collective pitch, and for M t =0.877 and collec- 
tive pitch angle 6 =8° in figure L 25. Both cases 
indicate good agreement; in the lifting case of 
figure 25, further improvement can be achieved by 
refining the wake model. 

For forward flight, computed pressure distribu- 
tions are compared with experimental data at a lo- 
cation near the blade tip in figure 26 for M t =0.8, 
advance ratio y=0.2, and zero collective pitch. 
Again, good agreement is demonstrated. 

The long-term objective of this project is to 
include viscous effects in the calculation of tran- 
sonic multibladed rotor flowfields. 


* Agarwal, R. K. and Deese, J. E. , "An Euler Solver 
for Calculating the Flowfield of a Helicopter 
Rotor in Hover and Forward Flight," McDonnell 
Research Laboratories, St. Louis, M0, 1987. 


7. Transonic Wing/Body Calculations Using the 
Slender-Layer Navier-Stokes Approximation (R. K. 

Agarwal and J. E. Deese, McDonnell Douglas Research 
Laboratori es , St. Louis, MO) 

With the increasing availability of supercom- 
puters, the ability to compute transonic flowfields 
has progressed to solution of Reynolds-averaged 
Navier-Stokes equations for aircraft configura- 
tions. A viscous flow code, designated MDSSL30, 
has been developed at McDonnell Douglas Research 
Laboratories for calculation of transonic wing/body 
flowfields . 

The code computes the turbulent flowfield by 
solving the three-dimensional, Reynolds-averaged, 
slender-layer approximation to the Navier-Stokes 
equations on body-conforming, curvilinear grids. 

In the slender-layer approximation to the Navier- 
Stokes equations, only streamwise diffusion terms 
are neglected; diffusion terms in the other direc- 
tions are retained. Therefore, slender-layer equa- 
tions are suitable for calculating spanwise separa- 
tion and the flow along a wing-body junction. 

These equations are solved by employing Jameson's 
finite-volume explicit Runge-Kutta time-stepping 
scheme. The calculations performed so far have 
used the Baldwin-Lomax turbulence model. 

The details of the methodology and the test 
calculations performed to verify the code are de- 
scribed by Agarwal, Deese, and Underwood (1985) and 
Agarwal and Deese (1984). Recently, fine-grid cal- 
culations for turbulent flow over an 0NERA-M6 wing 
were performed on a Cray X-MP/48; the code was 
fully optimized and mi crot asked for peak perfor- 
mance on the f our-processor Cray X-MP/48, as docu- 
mented by Booth and Misegades (1986). 
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Fig. 24 Surface pressure distributions on a nonlifting rotor in 

hover; M t = 0.52, 0 C = O°, AR = 6, NACA 0012 blade, and 
97x33x21 mesh. 


Computed results at various spanwise stations 
of the 0NERA-M6 wing are compared with experimental 
data in figure 27, with good agreement. The code 
has also been exercised by Deese and Agarwal* for 
calculation of fighter- and transport-type wing/ 
body flowfields; excellent global grids were gener- 
ated with the three-dimensional procedure of Chen, 
Vassberg, and Peavey (1985). Surface grid-lines of 
a transport wing/body configuration are shown in 
figure 28. Comparisons of calculated pressure dis- 
tributions at various spanwise stations with 
experimental data are shown in figure 29; the 
agreement is again good. 

The long-term objective of this project is to 
calculate turbulent flowfields for increasingly 
complete transport and fighter aircraft configura- 
tions. The project has also been accepted for ac- 
cess to the NAS Cray-2. 



Fig. 25 Surface pressure distributions on a lifting rotor in hover; 
M, = 0.877, 0 C = 8° , AR = 6.0, NACA 0012 blade, and 
97x33x21 mesh. 


8. Concluding Remarks 

An extensive program to develop advanced CFD 
codes is being conducted within the components of 
the McDonnell Douglas Corporation for applications 
to helicopters, transport and fighter aircraft, and 
missiles and hypersonic vehicles. Efficient use of 
large computers, including multiple-processor fa- 
cilities, is receiving special attention, and ac- 
cess to the NAS facility in several of these areas 
is greatly appreciated. 


*Deese, J. E. and Agarwal, R. K., "Navier-Stokes 
Calculations of Transonic Viscous Flow About 
Wing-Body Configurations,” McDonnell Douglas 
Research Laboratories, St. Louis, MO, 1987. 
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Fig. 29 Navier-Stokes solutions for flow past a transport- 
type wing-body; M a = 0.76, Re = 6.39 million, 

a = 2\ 
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